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Abstract. The main objective of this article is to study the dynamic tran- 
sition and pattern formation for chemotactic systems modeled by the Keller- 
Segel equations. We study chemotactic systems with either rich or moderated 
stimulant supplies. For the rich stimulant chemotactic system, we show that 
the chemotactic system always undergoes a Type-I or Type-II dynamic tran- 
sition from the homogeneous state to steady state solutions. The type of 
transition is dictated by the sign of a non dimensional parameter b. For the 
general Keller-Segel model where the stimulant is moderately supplied, the 
system can undergo a dynamic transition to either steady state patterns or 
spatiotemporal oscillations. From the pattern formation point of view, the 
formation and the mechanism of both the lamella and rectangular patterns 
are derived. 



1. Introduction 

Chemotaxis is a remarkable phenomenon occurring in many biological individ- 
uals, which involves mobility and aggregation of the species in two aspects: one 
is random walking, and the other is the chemically directed movement. For ex- 
ample, in the slime mould Dictyostelium discoideum, the single-cell amoebae move 
towards regions of relatively high concentration of a chemical called cylic-AAMP 
which is secreted by the amoebae themselves. Many experiments demonstrate that 
under some properly conditions a bacterial coloy can form a rather regular pattern, 
which is relatively stable in certain time scale. A series of experimental results on 
the patterns formed by the bacteria Escherichia coli (E. coli) and Salmonella ty- 
phimurium (S. Typhimurium) were derived in 2,3, where two types of experiments 
were conducted: one is in semi-solid medium, and the other is in liquid medium. 
They showed that when the bacteria are exposed to intermediates of TCA cycle, 
they can form various regular patterns, typically as ringlike and sunflowerlike for- 
mations. In all these experiments, the bacteria are known to secrete aspartate, a 
potent chemoattractant; also see [10} [T]. 

The most interesting work done by Budrene and Berg are the semi-solid experi- 
ments with E. coli and S. typhimurium. A high density bacteria were inoculated in 
a petri dish containing a uniform distribution of stimulant in the semi-solid medium, 
i.e. 0.24% soft agar in succinate. The stimulant provides main food source for the 
bacteria. In a few days, the bacteria spread out from the inoculum, eventually 
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covering the entire surface of the dish with a stationary pattern where the higher 
density population is separated by regions of near zero cell density. The S. ty- 
phimurium patterns are concentric rings and are either continuous or spotted; see 
Figure 1X7X1 The E. coli patterns are more complex with symmetry between individ- 
ual aggregates. A large number of patterns has been observed. The most typical 
forms are concentric rings, sunflower type spirals, radial stripes, radial spots and 
chevrons. In the process of pattern formation, the population of bacteria has gone 
through many generations. 




(a) (b) 

Figure 1.1. The black ring line and spot represent high density of bacteria. 



The liquid experiments with E. coli and S. typhimurium exhibit relatively simple 
patterns which appear quickly in a few minutes, and last about half an hour before 
disappearing. Two types of patterns are observed, and they rely on the initial 
conditions. The simplest patterns are produced when the liquid medium contains a 
uniform distribution of bacteria and a small amount of the TCA cycle intermediate. 
The bacteria collect in aggregates of about the same size over the entire surface of 
the liquid. The second type of patterns appears when a small amount of TCA is 
added locally to a special spot in a uniform distribution of bacteria. In this case, 
the bacteria begin to form aggregates which occur on a ring centered about the 
special spot, and in a random arrangement inside the ring. In particular, in these 
liquid experiments, the timescale to form patterns is less than the time required for 
bacterial birth and death. Therefore, the growth of bacteria does not contribute to 
the pattern formation process. 

Here we have to address that in these experiments, none of the chemicals placed 
in the petri dish is a chemo-attractant. Hence, the chemoattractants, which play 
a crucial role in bacterial chemotaxis, are produced and secreted by the bacteria 
themselves. 

In their pioneering work [5], E. F. Keller and L. A. Segel proposed a model in 
1970, called the Keller-Segel equations, to describe the chemotactic behaviors of 
the slime mould amoebae. In their equations, the growth rate of amoeba cells was 
ignored, i.e., the model can only depict the chemotaxis process in a small timescale, 
as exhibited in the liquid medium experiments with E. Coli and S. Typhimurium 
by [11 [3]. However, in the semi-solid medium experiments, the timescale of pattern 
formation process is long enough to accommodate many generations of bacteria. 
Therefore, various revised models were presented by many authors, taking into 
consideration the effects of the stimulant (i.e. food source) and the growth rate of 
population; see among others [TU] and the references therein. Also, there is a vast 
literature on the mathematical studies for the Keller-Segel model; see among others 

Elan nam]. 

The main objective of this article is to study the dynamic transition and pattern 
formation for chemotactic systems modeled by the Keller-Segel equations. The 
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study is based on the dynamic transition theory developed recently by the authors. 
The key philosophy for the dynamic transition theory is to search for all transition 
states. The stability and the basin of attraction of the transition states provide 
naturally the mechanism of pattern formation associated with chemotactic systems. 

Another important ingredient of the dynamic transition theory is the introduc- 
tion of a dynamic classification scheme of transitions, with which phase transitions 
are classified into three types: Type-I, Type-II and Type-Ill. In more mathemat- 
ically intuitive terms, they are called continuous, jump and mixed transitions re- 
spectively. Basically, as the control parameter passes the critical threshold, the 
transition states stay in a close neighborhood of the basic state for a Type-I transi- 
tion, are outside of a neighborhood of the basic state for a Type-II (jump) transition. 
For the Type-Ill transition, a neighborhood is divided into two open regions with 
a Type-I transition in one region, and a Type-II transition in the other region. 

Two types of Kcller-Segel models are addressed in this article. The first is the 
model for rich stimulant chemotactic systems (with rich nutrient supplies) . In this 
case, the equations are a two-component system, describing the evolution of the 
population density of biological individuals and the chemoattractant concentration. 
In this case we show that the chemotactic system always undergoes a Type-I or 
Type-II dynamic transition from the homogeneous state to steady state solutions. 
The type of transition is dictated by the sign of a nondimensional parameter b. For 
example, in a a non-growth system in a narrow domain, for the spatial scale smaller 
than a critical number, the system undergoes a Type-I (continuous) transition, 
otherwise the system undergoes a Type-II (jump) transition, leading to a more 
complex pattern away from the basic homogeneous state. 

The second is a more general Keller-Segel model where the stimulant is mod- 
erately supplied. In this case, the model is a three-component system describing 
the evolution of the population density of biological individuals, the chemoattrac- 
tant concentration, and the stimulant concentration. In this case, the system can 
undergo a dynamic transition to either steady state patterns or spatiotemporal os- 
cillation. In both cases, the transition can be either a Type-I or Type-II dictated 
respectively by two nondimensional parameter 6 and b\. 

For simplicity, we consider in this article only the case where the first eigenvalue 
of the linearized problem around the homogeneous pattern is simple (real or com- 
plex), and we shall explore more general case elsewhere. In the case considered, 
for the Type-I transition, when the linearized eigenvalue is simple, we show that 
both the lamella and rectangular pattern can form depending on the geometry of 
the spatial domain. Namely, for narrow domains, the lamella pattern forms, other- 
wise the rectangular pattern occurs. For Of course, for Type-II transitions, more 
complex patterns emerge far from the basic homogeneous state. 

The paper is arranged as follows. Section 2 introduces the Keller segel model. 
The rich stimulant case is addressed in Section 3, and the general there-component 
system is studied in Section 4. Section 5 explores some biological conclusions of the 
main theorems. 



4 



MA AND WANG 



2. Keller-Segel Model 
The general form of the revised Keller-Segel model is given by 



(2.1) 



dui . 
— - = fciAui - xv(ui Vh 2 + aim 
at 

du 2 

~dt 

du 3 

dt 



Q'O + U 3 



k 2 Au 2 + rim - r 2 u 2 , 
k 3 Au 3 - r 3 mu 3 + q(x), 



where m is the population density of biological individuals, u 2 is the chemoattrac- 
tant concentration, u 3 is the stimulant concentration, q(x) is the nutrient source, 
and x is a chemotactic response coefficient. 

Equations (|2.1|) arc supplemented with the Neumann condition: 



(2.2) 



d{ui,u 2 ,u 3 ) 
dn 



on 90. 



For simplicity, we consider in this article the case where the spatial domain is 
a two-dimensional (2D) rectangle: 

0= (Q,h) x (0,/ a ) for h^l 2 . 

It is convenient to introduce the nondimensional form of the model. For this pur- 
pose, let 



(2.3) 



t = t'/r 2 , x = \Jk 2 jr 2 x , 

ui = < l /a~2~u' 1 , u 2 = k 2 u' 2 /x, u 3 =a u 3 , 



and we define the following non dimensional parameters: 

\ = riy/a2~x/ r 2k 2 , a = a 1 a 2 /r 2l (i = ki/k 2 , 
r = k 3 /k 2 , S = r 3y /a 2 ~/r 2 , 5 = q/r 2 a . 



(2.4) 



Then suppressing the primes, the non-dimensional form of the Keller-Segel model 
is given by: 



(2.5) 



dm . . 

— — = fi Am — V iui vu 2 ) + am 
at 

du 2 

~dt 

du 3 

dt 
du 
On 



"3 



1 + M3 



1 ' 



Au 2 - u 2 + Ami, 
rAu 3 - 5mu 3 + 5 , 
= 0, 



an 

u(0) = uq in Q. 
The non-dimensional of Q is written as 

Q = (0, L x ) x (0, L 2 ) with Li ^ L 2 . 
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Often times, the following form of the Keller-Segel equations is discussed in some 
literatures: 



(2.6) 



dui , . 

at 
du 3 

at 

- Au 2 + u 2 — Xui. 
du 



«3 



1 + M3 



1 ' 



rAu 3 - Suiu 3 + S , 



dn 



= 0, 



an 

u(0) = u . 

The biological significance of (|2.6I) is that the diffusion and degradation of the 
chemoattractant secreted by the bacteria themselves are almost balanced by their 
production. The main advantage of (|2 .6[) lies in its mathematical simplicity, and 
as we shall see from the main results of this article, the the main characteristics of 
the pattern formation associated with the model are retained. 



3. Dynamic Transitions for Rich Stimulant System 



3.1. The model. We know that as nutrient u 3 is richly supplied, the Keller-Segel 
model (|2.1[) is reduced to a two-component system: 

dm 

~dt 

du 2 

(3.1) ~dt 

d{ui,u 2 ) 



dn 
u(0) = u Q 



= /iAui — V(uiVu2) + aui(l — u x ), 
= Au 2 — u 2 + Xui, 
= 0, 



on 



It is easy to see that u* = (1, A) is a steady state of ()3.1|) . Consider the deviation 
from u*: 

u — u + u . 

Suppressing the primes, the system p.ip is then transformed into 



(3.2) 



— — = fiAux — 2au\ — Au 2 — V(uiVu2) 
at 

— - = Au 2 -u 2 + Xui, 
at 



1) 



d(u 1 ,u 2 ) 



dn 
u(0) = M 



0. 



an 



3.2. Dynamic transition and pattern formation for the diffusion and degra- 
dation balanced case. We start with an important case where the diffusion and 
degradation of the chemoattractant secreted by the bacteria themselves are almost 
balanced by their production. In this case, the second equation in (|3.2|) is given by 



= Au 2 — u 2 + Atti. 



6 



MA AND WANG 



With the Newman boundary condition for U2, we have U2 = [—A + 1] ui and the 
functional form of the resulting equations are given by 

(3.3) <^± = Cx u 1 + G(u 1 ,\), 

ot 

where the operators C\ : Hi — > H and G : H\ x M — > K are defined by 

C\U\ = /iAui — 2aui — AA[— A + I\~ ui, 
l ' > ' 4 ' G(wi,A) = -AV(uiV[-A + J]-V)-3cm 2 -mt?. 

Here the two Hilbert spaces H and H\ are defined by 

H = L 2 (fl), H l = {u x e H 2 (fl) I = on fi}. 

an 

To study the dynamic transition of this problem, we need to consider the lin- 
earized eigenvalue problem of (|3.3|) : 

(3.5) C x e = p{\)e. 

Let pk and e k be the eigenvalues and eigenfunctions of —A with the Neumann 
boundary condition given by 

kiirxi k^isx-i 2 ( k% 



(3.6) e fe =cos— cos— , p k = tt 72+72 

L\ L2 ^2 

for any k = (k±,k2) € N 2 ,. Here N+ is the set of all nonnegative integers. In 
particular, eo = 1 and p = 0. 

Obviously, the functions in (|3.6j) are also eigenvectors of (|3.5j) . and the corre- 
sponding eigenvalues /3fc are 

(3.7) /? fe (A) = - m -2a+-^. 
Define a critical parameter by 

/ Q Q \ \ • (gfc + jOQfgfc + 20! ) 

(3.8) A c = mm . 

Let 

S = {K = [Kx^K^) G achieves the minimization in (|3.8[) } . 
Then it follows from (JHIZJ) and that 
< if A < A c 

(3.9) Pk(\){ =0 if A = A c VK= (K 1 ,K 2 ) eS, 
> if A > Ac 

(3.10) /3 fc (A c )<0 Vfc € Z 2 with fc ^ 5. 
Notice that for any X = {K 1 ,K 2 ) *E S, K =/= 0, and 

(3.11) Ae= (px + l)(^ + 2a) 

PK 

We note that for properly choosing spatial geometry, we have 



(3.12) pK = 7T ^-± + -£j = A /— - VX=(if 1 ,X 2 )e«S, 

(3.13) A c = 2a + /i + 2 v / 2a/I. 
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Conditions (|3.9[) and (|3.10[) give rise to a dynamic transition of Q3.3P from (u, A) 
(0, A c ). For simplicity, we denote 

Ki±(Ki,0), K 2 ±(0,K 2 ), 

and K = {K \,K 2 ) € S. Also, we introduce a parameter as 

(<,^A^ h q ^ No 24-15sign(^ 1 X 2 ) - 

(3.14) 6 = - 3 /W + ^12 - 4 „ 2sign( ^ X2) 

(2fJfK + a)(2pp 2 K + 28ap g + 4a - /xgg) 

[1 + sign(X 1 X 2 )] ■ [(jjipaK + 2a)(l + /j 2 k) - A c p 2Jf ] 

2(2 fip K p Kl + 40:/?^ - 3ap_R-) 

(1 + signal )p 2 K [(np2K 1 + 2a)(l + p 2 ifj - A c p 2 / fl ] 
x [(/ipif + 2a)(2p|- 1 - &p Kl PK 2 - Pk) + 6ap K {4:pK 1 + 1)] 

2(2pp K pK 2 + 4apK 2 - 3apjf ) 

(1 + signK 2 )p 2 K ((pp2K 2 + 2a)(l + p 2i f 2 ) - \ c p 2 K 2 ) 
x [(/i/Jjf + 2a)(2p|- 2 - 6p Kl PK 2 - px) + 6ap K (4:p K2 + 1)]. 

The following is the main dynamic transition theorem, providing a precise crite- 
rion for the transition type and the pattern formation mechanism of the system. 

Theorem 3.1. Let b be the parameter defined by \3. 1J$ . Assume that the eigen- 
value Pk satisfying i3. 9\) is simple. Then, for the system VS. 3]) we have the following 
assertions: 

(1) The system always undergoes a dynamic transition at (u, A) = (0, A c ). 
Namely, the basic state u = is asymptotically stable for A < X c , and 
is unstable for A > A c . 

(2) For the case where b < 0, this transition is continuous (Type-I). I particular, 
the system bifurcates from (0, A c ) to two steady state solutions on A > A c , 
which can be expressed as 



±, u ,1 PkW K 1 ttx 1 K 2 7TX 2 f a l/2\ 
(3.15) ui (x, A) = ±- ] j — cos _ cos _ + j , 

and u 1 (x, A) are attractors. 
(3) For the case b > 0, this transition is jump (Type-II), and the system has 
two saddle-node bifurcation solutions at some A*(0 < A* < A c ) such that 
there are two branches and u 2 of steady states bifurcated from (v* , A* ) , 
and there are two other branches and bifurcated from (u* , A* ) . In 
addition, v* and v£ are saddles, i> 2 and v% are attractors, with ,v$ 
as A — > A c . 

Two remarks are now in order. 

Remark 3.1. From the pattern formation point of view, for the Type-I transition, 
the patterns described by the transition solutions given in (|3.15[) are either lamella 
or rectangular: 

lamella pattern for K\K 2 = 0, 

rectangular pattern for K\K 2 7^ 0. 

In the case where b > 0, the system undergoes a more drastic change. As A* < 
A < A c , the homogeneous state, the new patterns v 2 and v± are metastable. For 
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A > A c , the system undergoes transitions to more complex patterns away from the 
basic homogeneous sate form. 

Remark 3.2. If we take the growth term f(u) as / = au\{l — U\) instead of 
/ = au\{l — u\) in p.ip . p. 21) and p.3p , then Theorem 13. II still holds true except 
the assertion on the existence of the two saddle-node bifurcation solutions, and the 
parameter should be replaced by 

{2pp K + a)(2(i\ 2 K + 10ap K + a - fxp K ) 



pp K + a 



2(1 + signK 1 K 2 )[(pp2K + a)(l + p 2 x) - Kp2K\ 
(2pp K p Kl + 2ap Kl - ap K ) 



(1 + signK 1 )X 2 K [(pp2K 1 + ac)(l + p 2 K 1 ) - Kp2K x \ 
x [{Wk + a)(2X% 1 - 6p Kl PK 2 - Pk) + 2ap K (Ap Kl +1)] 
(2p,pKPK 2 + 2apK 2 - apK) 

(1 + SlgnK2)\ 2 K [{pp2K 2 + Ct)(l + P2Kt) - X c p2K 2 \ 

x [{ppK + a)(2\ 2 Ki - 6p Kl PK 2 - Pk) + 2ap K (Ap K . 2 + 1)]. 

3.3. Pattern formation and dynamic transition for the general case. Con- 
sider the general case (|3.ip . In this case, the unknown variable becomes u = (ui, M2), 
and the basic function spaces are then defined by 

(jit 

H = L 2 (n,R 2 ), H x = {u E H 2 (n,R 2 ) I — = on fi}. 

on 

Let L\ : Hi — > H and G : Hi -> H be defined by 
(3.16) 



, , aA-2a -A 1 
L\u = [ x A _ x I u, 



n i \ _ ( — V(miVu2) — 3au 2 — au\ 

Lr(U) — I q 

The linearized eigenvalue problem of (|3.2|l is 
(3.17) L xV = p Vl 

where L\ : Hi — > H is defined by (|3. 16|) . Let be the matrices given by 

(3.18, S i = (-C« + 2<1) -(„f + l) 

where p k are the eigenvalues as in p. 61) . It is easy to see that all eigenvectors ifk 
and eigenvalues Pk of (|3.17p can be expressed as follows 



(3.19) ip k 



Cfei e fc 

£k2^k 



(3.20) Bl ( l kl ) = (3 k ( i; 1 



Cfe2 / V £ 



s/,'2 



where e k are as in (|3.6p , and Pk ar e also the eigenvalues of B^. By p.!8p . p k can 
be expressed by 



PkW = \\- B± \/S 2 - 4((p fe + l)(pp k + 2a) - \p k ) 



(3.21) ^ fe w 2 

B = (/i + l)p fe + 2a + l. 
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Let A c be the parameter as defined by (13. 8[) . It follows from ()3.2ip and (|3.8p that 

f <0 if A < Ac, 
(3.22) Pk( X ){ = if A = A c , 

[ > if A > Ac, 

r o 9 ^ / Re/\7(A C )<0 VfceZ 2 , 

[ ' \Re£+(A c )<0 Vfc e Z 2 with Pk ^ PK , 

with K = {Ki,K 2 ) as in (|3~TT|) . 

Then we have the following dynamic transition theorem. 

Theorem 3.2. Let 6 be the parameter defined by \3. 11$ . Assume that the eigen- 
value satisfying 113. 22\) is simple. Then the assertions of Theorem \3.1\ hold true 
for 13. ty) . with the expression 113.15]) replaced by 

u x =±y'a/3+(A)l p ^ I cos— — cos— — + o{\/3+\ ' ), 

8(/ip g + + 2a + 1) 
( P K + 1) 3 |&| 

3.4. Proof of Main Theorems. 



Proof of Theorem \3.1\ Assertion (1) follows directly from the general dynamic tran- 
sition theorem in Chapter 2 of J7]. To prove Assertions (2) and (3), we need to 
reduce (|3. 31) to the center manifold near A = A c . We note that although the under- 
lying system is now quasilinear in this general case, the center manifold reduction 
holds true as well; see [6] for details. 

To this end, let u = set + $, where $(x) the center manifold function of (|3.3[) . 
Since L\ : H\ — > H is symmetric, the reduced equation is given by 



(3.24) 

where G : Hi 
(3.25) 



dx 



Pk{\)x 



1 



(ex, ex) 
H is defined by (pT4|) . and 

(ex, ex) 



(G{xe K + $,\),e K ), 



e%dx= 2 - S ^ K ^M 



It is known that the center manifold function satisfies that $(x) = 0(x 2 ). A 
direct computation shows that 

(3.26) <G(xe K + <5>,\ c ),e K > 



= —ax 



e K dx — 6ax 



e K $dx 



+X c x / [e K Ve k ■ V(-A + + $Ve fe ■ V(-A + I)~ l e K ]dx + o(x 3 



It is clear that 



(-A+/rv 



PK + l 



ex, 



Ae K = -Pk&k- 
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We infer from ([3^6]) that 
(3.27) (G(xe K + <Z>,\c),e K ) 



e%dx — 6ax 



/ e 2 K $dx 
Jq 



+ X c x 



_PK + 1 

+PKe 2 K (-A + I)' 1 ®] dx + o(x 3 ). 

Using the approximation formula for center manifold functions given in (A. 11) 
in [5], $ satisfies the equation 



(3.28) -L Ac $ = G 2 (xe K ,Xc)+o(x 2 ) 

PK + 1 



3a I e K — 



A, 



PK + 1 



+ o(x 2 ). 



In view of (13. 61) . we find 



(3.29) 



e 2 K = t[ 6 o + e 2K\ + e 2 K 2 + em], 



|Ve_R"| 2 = -[pKeo + {pK 2 - PK^)Z1K X + {PK X - PK 2 )e 2 K 2 ~ Pk&2k\- 



Thus, l|3.28[) is written as 
(3.30) -L Ac $ 



x 

"T 



/ 2p Kl X c 
—6ae + \ — — oa e 2 K 1 



+ 



2pK 2 X c 

PK + 1 



\PK + l 

- 3a I e 2 K 2 



\PK + l 



e 2K ~ 3a e 2K 



+ o(x 2 ). 



Denote by 
(3.31) 
Note that 
(3.32) 



$ = $ e + & 2 K 1 e 2 K 1 + $2K 2 e 2 K 2 + ^2K^2K- 



L\ c e 2K 



1 



1 + P2K 

Then, by l|3~TT|) and (1330)) - (13321) we obtain 

3 



[(1 + P2k){pPk + 2a) - \cP2k] e 2K - 



(3.33) 



$o = - 

®2K 1 = 
<S>2K 2 = 
$>2K = 



(1 + P2K 1 )(2fJ.pKPK 1 + ^apK 1 ~ 3apg) 
4pif[(l + P2K 1 ){PP2K 1 + 2a) - P2K 1 K] ' 

(1 + P2K 2 ){2pp K p K2 + 4ap K2 - 3ap g ) 

4pif[(l + P2K 2 )(PP2K 2 + 2a) - P2K 2 X C ] ' 

(1 + P2k)Qwk + a) 
4[(1 + P2k)(pP2K + 2a) - p 2 A'A c ] ' 
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Inserting (|3~3Tj) and fl3l| into ([3^7]) we get 

(3.34) < G[xe K + $, A c ), e K >= -ax 3 / e^cfa; 

Jo 

6qx(2 - signal if 2 )) /" r , 2 , A 2 lfh 2 , * 21. 
4 J [®o e o + $2x^2^ + ®2K 2 e 2K2 + $2Ke 2K \ dx 

[pK^oel + (pk 2 - PK 1 )<S>2K 1 e 2 2Kl 



\ c x{2 - sign(X 1 ^ 2 )) 
4(p K + 1) 

A c :r(2 - sign(X 1 ^ 2 )) 



^ 2 , PA 2 ~ PKl = 2 

PA$oe + — — ; — ; < ^2K 1 e 2Kl 



PKi ~ PK 2 ^ 2 

+ -T— $2K 2 e 2K2 

1 + P2if 2 

\ c p K x(2 - signal if 2 )) 



1 + P2Ki 

-3 > 2Are 2K 



1 + P2iC 



ri.r 



4> 



2Ki 



2K 2 2 
( 9 



t : c 2K -r — — — c 2K 

1 + P2K1 t + p 2 ^ 2 



2/\ 



1 + P2K 



'2K 



-ax 



e A K dx 



dx + o(x 3 ) 
fl\x(2 - sign^i^)) 



+ 



1 

1 + signif 2 
1 



1 



1 + signal 



A c (Pa 2 - PKi) , 1\ c p Kl 



I + PK 



1 + P2Kx 



6a §2Kt 



Ac(PA-x - PK 2 ) . 2X c pK 2 



1 + PK 1 + p2Ki 

KPK , 2\ c p K 



- 6a <&2K 2 



2(1 + sign(X 1 X 2 )) V 1 + PK 1 + P2K 
Also, we note that 



6a &2K 



dx + o(x 3 ). 



Li 



et-dxi / ej io dx 2 



24- 15sign(iVi.Fs: 2 ) 



G4 



Then, putting ([3~33"|) into ([3~34"|) we get 



(3.35) 



< G{xe K + $,A c ),e K >= 



(2-sign(^ 1 K 2 ))|r>|x 3 
32 



b + o(x 3 ), 



where b is the parameter given by (|3.14p . 

By (|3.24[) and p.35[) . we derive the reduced equation on the center manifold as 
follows: 



(3.36) 



^ = /3 K (X)x + -x 3 + o(x 3 ). 



Based on the dynamic transition theory developed in Chapter 2 in [7], we obtain 
Assertions (2) and (3), except that two saddle-node bifurcations occur at the same 
point A — A*. To prove this conclusion, we note that if u*(x) is a steady state 
solution of (13.31). then 



V*(x) = U*{x + 71") — U*(x — 7r) 

is also a steady state solution of (|3.3|) . This is because the eigenvectors (|3.6|) form 
an orthogonal base of H\. Hence, two saddle-node bifurcations on A < A c imply 
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(3.37) — = [3+(X)x + — - - < G(x ■ ifiK + <Pk >, 



that they must occur at the same point A = A* . Thus the proof of the theorem is 
complete. □ 

Proof of Theorem\EM Assertion (1) follows from ([3~2"2"j) and ([3~23")1 . To prove As- 
sertions (2) and (3), we need to get the reduced equation of (|3.2[) to the center 
manifold near A = A c . 

Let u = x ■ (fK + where ifiK is the eigenvector of (|3.17[) corresponding to (3k at 
A = A c , and &(x) the center manifold function of Q3.2p . Then the reduced equation 
of {HU read 

dt , KK ' <tp K ,tp* K > 

Here ip* K is the conjugate eigenvector of (px. 
By (13.191) and (|3.20p . <p K is written as 

(3.38) (p K = (fie*, ii&Kf \ 

with (£i,£ 2 ) satisfying 

(" ( " p r 2o) >=»• 

from which we get 

(3-40) (6,£a) = 0*- + l, Ac). 

Likewise, y|f is 

(3-41) ¥>x = (£e fc , ^e fc ) T , 

with (£* , £2 ) satisfying 

-(fxp K +2a) A c \( Ci 

which yields 

(3-42) (&&) = (pk+i,Pk). 

By (13.161) . the nonlinear operator G is 

G(ui,u 2 ) = G 2 (ui,u 2 ) + G 3 (ui,u 2 ), 



G 2 (ui, u 2 ) = — (VwiV«2 + uiAm 2 + 3a«i) 
G 3 (ui, u 2 ) = ~au 3 



1 


It is known that the center manifold function 

$(a;) = ($1(1), $ 2 (a;)) = 0(a; 2 ). 
Then, in view of (|3.38p and (|3.40[) . by direct computation we derive that 
(3.43) (G(x^e K + ^i, x^e K + *a), 

= / K2*i|VeA-| 2 -^iA* 2 4:-6aCi*ie|-]da; 

-aQ£,fx 3 / e^ete + o(x 3 ). 
Jn 
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Using the approximation formula for center manifold functions given in (A. 11) 
in [8], $ = ($1,^2) satisfies 



(3.44) -L Ac $ = -x 2 G 2 (^e k ,^e K ) + o(x 2 ) 

= -x 2 (a6|VeK| 2 + (3a£ 2 - te2 PK )e%) ( J 

From (13.61) we see that 







e2Ki)(l + e-2K 2 ) = 


1 , 


+ e 2 ic), 


e 2 K 










- e 2K t ){^ + e2K 2 ) 


+ (1 + 62,0(1 


- e2^ 2 ) 




PK 


. Pk 2 - PKi 


PKi - PK 2 


PK 




~e - 


f 4 e 2Kl - 


I" 4 e 2/<2 


—e 2K . 



Thus, (|3.44l) is written as 
fix 2 

(3.45) L Ac $ = - ^L-(3c«ieo + (3o& - %pjr 1 )e 2 jf a 

+ (3a£i - 2£_ 2 pK 2 )e2K 2 + (3a£i - 2£_ 2 p K )e2K) ( q ^ + o(x 3 ). 

Let 

(3 - 46) ( £ ) = ( «° ) e ° + ( S ) e ^ + ( S ) e2A " 2+ ( ¥ ) e ^ 

It is clear that 



*§ J efc ~ Bfc I *s 

where is the matrix given by (|3.18p . Then by (|3.45[) and (|3.46p we have 



$f \ = (3<-2a$2PfcK i? -i ( 1 



for k = K,K 1 ,K 2 , and B 2fc = . 
Direct computation shows that 



(3.47) 
(3.48) 
(3.49) 
(3.50) 



$0 \ _ 3£ 2 :r 2 / 1 
>°2 



y 8 v a c 



$ 2 ^i \ = g 1 (3ag 1 -2g 2PjCl ) / 1 + P2K1 
$ 2if2 y 4detB 2 /f 1 V A c 

<$>\ K * \ = Cl(3<l-26PK 2 ) / 1 + P 2K 2 

§l K2 ) 4det£? 2 /f 2 v A c 

<$>\ K \ = ^i(3agi-26px) / 1 + p 2 x 
$i A ' 7 4dctB2K I A c 
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Inserting ([OTjl into (p£33jl . by (HQO) and (pT42|) we get 

(2-signif 1 if 2 )(p x + l)|0| 



8a(pK + l) 3 x 3 J n e K dx 



{2-sigaK 1 K 2 )\n\ 
2 

+ l + signK 1 (UPK2 ~ PKl) " ^O*?* 1 * 

+ l + signg 2 (6(pgl " PK > ] ~ 6a ^ 2x 

2 - signify , c L _ t u2K 
2 + 6a§i)*i a; 

1 + signiv ! 1 + signif 2 



(p K + l)/?2A-(2 - signal K 2 ) , 2 a: 



+ o(x 3 ). 



By definition, we have 

PK X + PK 2 = PK, PlK = 4p K VK = (Ki,K 2 ), 

< ip, ip* >= [{p K + l) 2 + Pa- A c ] / e 2 K dx 

Jn 

= 2 _^lM {pK + 1)(/ipK + pK + 2a+ 1)|0|. 



(3.51) 



In view of (|3.47|) - (|3.50p . the reduced equation (|3.37[) is given by 
dx (pk + l) 3 bx 3 



dt HKy ' 8(np K + p K + 2a + 1) 



+ o(x 3 ), 



where b is the parameter as in (|3.14[) . Then the theorem follows readily from (|3.51|) . 
The proof is complete. □ 



4. Transition of Three-Component Systems 

4.1. The model. Hereafter 5q > is always assumed to be a constant. Hence, 
(|2.6[) has a positive constant steady state u* given by 



(4.1) («i,u5,«S) withuj = l- 



1/2 



u 2 =Xu 1 , u 3 u 1 = —. 



It is easy to see that is the unique positive real root of the cubic equation 



/ <h\ f So 



Consider the translation 

(4.2) {ui,u 2 ,u 3 ) -t (u{ + Ui,u 2 +u 2 ,u{ + ui). 
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Then equations (|2.6|) are equivalent to 

dm . „ 2 „ . aul , . 

— = /iA?ti - 2au x ui ~ u 1 Au 2 + ^ + u *^2 Ua + 9W, 

du 

-^j- — rA«3 — Su*v,3 — Su^ui — SU1U3, 

(4.3) 

— AU2 + U2 = Aui, 
d(ux,u 2 ,u 3 



= 0, 

an 



On 
u(0) = u Q , 

where u = (u l7 w 3 ), u 2 = A[— A + 1] — U\, and 
(4.4) g(u) = -V(uiVu 2 ) - iau*iu\ - au\ + 
au*u 3 au\u 3 au^ui 



g(ui + up(u 3 + m|) 

1 + U3 + M 3 



l + u* 3 (I + U3) 2 I + U3' 
The Taylor expansion of g at u = is expressed by 

/ \ , , „ 2 a Ul u 3 au\u% 

g(u) = -V^-S^i^-^ 

Let 

"2/o 1D>2\ 



i?i = {w G 

Define the operators L\ : Hi — > H and Ga : i?i — > H by 
(4.5) 



£Ti = {m G H 2 (Sl,R 2 )\ -H = on 90}. 

on 



L AU = [ ^ A " 2aU * 2 " M A[-A + J]- 1 \ / Ul 

I rA - <5u| /V u 3 



G(u,X) = ( 9 , {U) ) 



Then the problem (14.31) takes the following the abstract form: 

du 

(4.6) Tt =L x u + G{uAl 

u(Q) = Mo- 
lt is known that the inverse mapping 

[-A + /]- 1 :H^Hi 
is a bounded linear operator. Therefore we have 

L\ : H\ — > H is a sector operator, and 

Ga : He -> i? is C°° bounded operator for # > -. 

We note that the transition of (|4.3|) from it = is equivalent to that of (|2.6|) from 
u = u . 
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Theorems 13.11 and 13.21 show that a two-component system undergoes only a 
dynamic transition to steady states. As we shall see, the transition for the three- 
component system l|2.5p is quite different - it can undergo both steady state and 
spatiotemporal transitions. 

4.2. Linearized eigenvalue of (|2.6p . The eigenvalue equations of (|2.6p at the 

steady state (u^it^Mg) given by (|4.1[) in their abstract form are given by 

(4.7) L xV > = pip, 

where L\ : H\ — > H as defined in (|4.5|) . The explicit form of (|4.7j) is given by 

/ pA-2auf -A^AI-A + J]- 1 jgfa, 
^ -5u* 3 rA- 3 5u\ 

As before, let pk and ek be the eigenvalue and eigenvector of —A with Neumann 
boundary condition given by (|3.6|) . and let 

Then, it is easy to see that tp k is an eigenvector of (|4.7p provided that (£fci> £fc3) <= ^ 2 




is an eigenvector of the matrix Ai' 



At(f 1 )=Pk(f 1 
with 



(4.8) A 



fc 



— <5«3 — rpfc — 5u\ J 

The eigenvalues of A^ , which are also eigenvalues of (|4.7[) , are expressed by 

1 



J k 

(4.9) 



o± a /a 2 -4detA£ 



,\ ^Pk U l , ; > 

trA fc = — ^p/c - 2cmi - rp fc - 0%. 

1 + Pfc 



To derive the PES, we introduce two parameters as follows: 
(4.10) A c = min i£E±A \ Wk + 2a uf + rp K + 6u*,] , 

PK p K U\ 



(4.11) Ac = min^±±l 

PK PKU X 



i 1 *2 , a(5 

p,pK + 2au x 1 



(l + u|) 2 (rp K +<K). 
Let K — (KijKs) and K* = (K^K^) be the integer pairs such that px and pk* 
satisfy (|4.10l) and (I4.11j) respectively. 

Theorem 4.1. Lei A c and A c 6e i/ie parameters defined by and re- 

spectively. Then we have the following assertions: 

(1) As A c < A c , the eigenvalues P K (X) of J^.ffp are a pair of conjugate complex 
numbers near A = A c , and all eigenvalues of J^.ff[ ) satisfy 

( <0 A<A C , 

(4.12) i?e/?|(AW =0 A = A C , 

I >0 A>A C , 

(4.13) Re p£(A e ) < 0, Vfc € Z 2 mift p fe ^ p A - 
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(2) As A c < A c , the eigenvalue /3^«(A) is real near A = A c , and all of &4-9\ ) 
satisfy 

< 0, A < A c 

(4.14) P K .W{ =°» X = X < 



(4-15) 



> 0, A > A c , 
Re/3+{\ c ) < 0, Vfc G I? with p k ^ p K , 
RePkW)<0, V|fc|>0. 



Proof. By (|4.9p we can see that /3^(A) are a pair of complex eigenvalues of (j4.7|) 
near some A = A*, and satisfy 

f <0, A < A*, 
Re/3±(AW =0, A = A*, 
[>0, A > A*, 

if and only if 

tr^lf = 0, det Al* > 0. 
Likewise, /?^(A) is real near A = A* and satisfies 

f <0, A < A*, 
=0, A = A*, 
[ > 0, A > A*, 

if and only if 

trAf < 0, det A£ = 
Due to the definition of A c and A c , when A c < A c we have 



(4.16) 



and when A c < A c , 



tr A K " = 0, 

tr A% c < 0, Vfc G Z 2 with p k ± p K , 
det Al" > 0, V|fc| > 0, 



det Ab. = 0, 



(4.17) det Al" > 0, Vfc G Z 2 with pfc ^ p A ', 

tr Al" < 0, V|fc| > 0. 

It is known that the real parts of fit are negative at A if and only if 

det Al > 0, tr Al < 0. 

Hence, Assertions (1) and (2) follow from (|4.16l) and (14. 17)) respectively. The the- 
orem is proved. □ 

4.3. Dynamic transition theorem for (|2.6[) . Based on Theorem 14. 11 we imme- 
diately get the following transition theorem for (|4.3I) . 

Theorem 4.2. Let A c and A c be given by and J^. _?_?[ ) respectively. Then the 

following assertions hold true for OTM); 
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(1) When A c < X c , the system undergoes a dynamic transition to periodic so- 
lutions at (u, A) = (0,A C ). In particular, if the eigenvalues fi K satisfying 
{4-.12\j are complex simple, then there is a parameter bo such that the dy- 
namic transition is continuous (Type-I) as bo < 0, and is jump (Type-II) as 
bo > with a singularity separation of periodic solutions at some A* < A c . 

(2) When A c < A c , the system undergoes a dynamic transition to steady states 
at (u,X) — (0,A C ). ///3^»(A) satisfying H4-14\ ) is simple, then there exists a 
parameter b± such that the transition is continuous as b\ < 0, and jumping 
as b\ > with two saddle-node bifurcations at some A < A c from (u + ,X) 
and (u~, A). 

Remark 4.1. By applying the standard procedure used in the preceding sections, 
we can derive explicit formulas for the two parameters bo and b\ in Theorem 14.21 
However, due to their complexity, we omit the details. Instead in the following, we 
shall give a method to calculate bo, and for b\ we refer the interested readers to the 
proof of Theorem 13.21 

4.4. Computational procedure of bo- The procedure to compute the parameter 
bo in Assertion (1) of Theorem 14.21 is divided into a few steps as follows. 



Step 1. The reduced equations of (|4.6[) to center manifold at A = A c are 
expressed by 

dx 1 
(4.18) * «*" > 

at < ip, ip* > 

where ip and ip are the eigenvectors of L\ at A = A c , ip* and ip* the conjugate 
eigenvectors, and L\,G\ : Hi — > H the operators defined by (|4.5p . <& is the center 
manifold function. 

Step 2. Solving the eigenvectors ip,ip and their conjugates ip*,ip*. We know 
that tpi and ip* are 

<P= (^iejr,6eK), ip = (vi^k, 172 ejc), 

and & , £* satisfy 

\ Su* 3 -r PK -5u\ / V £2 J ' \V2 

\ -5u* 3 -r PK -5u$ J \V2 J P V & 



and 



^f-PPK~2auf -Su* 3 \fCi\__Jvl 



(\ il n 



(i+«*) 

A c pk^i _ „„„ _ 9 ^„*2 
1+PK 

(1+«|) 2 




'/2 
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where A c is as in (|4.10[) . efc as in (|3.6|) . and 



A c pjftti r, *2 r * 

- PPk — 2au 1 =rp K + ou 1 , 



l + PK 

Here, we use that ulul — So/6. From these equations we obtain 



p = det A% = f , ° ° - icfPK + 5u\f. 



£1 = -( r PK + Su*), 6 = Sus, 

(4.20) g = g = °' 

v ' q = o, e 2 - -p, 

T)l = 6u* 3 , nliPK + Su\. 
Due to P~lT)f and (g^ZD]) we see that 



< ¥,<f* >=< i'* >= (mVi +V21I2) / e K dx = -Spul / e K dx, 

Jn Jn 

< (p,ip* >=< ip,tp* >= 0. 

Step 3. We need to calculate 

< G(xip + yip + $, A c ),w* >, with = ^*,w 2 = V '■ 
By (143)1 we have G = G 2 + G 3 , and 

-1, ,.\ Q „,,.,*, .2 | Oh>lU)2 



G 3 (w,A) = I 1 (!+«3) 3 (i+«3) 4 J , 

for w = (wijWa) € H 1 , By (|4.20l) we find 

< G 3 (xtp + + A c ), <p* >= 0. 

Noting that 



we have 



exejeidx = 0, VAT, J, / S Z 2 , 
$ = ($ 1 ,$ 2 ) = 0( a; 2 ), 

< G 2 (x^ + ?/?A + $,A c ),^*) 

/ [£*ejcff2i + Q, e K922]dx 
Jn 

(by £ = 0) 

/ Q e K[-5{x^ie K + yi]i e K + <&l){x& + OT2 + ®2)]dx 
Jn 

-8Q ( &x / <i>ie 2 K dx + £ix / $%e 2 K dx + / $ 2 e 2 K dx 
V is; > 
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Thus, we get 

(4.21) <G(xv + yip + $,A c ),<p* > 



= -6& 



&x / $ie 2 K dx + £iz / <i> 2 e 2 K dx + my / $ 2 e 2 K dx 
Jn Jn Jn 



o(3). 



In the same fashion, we derive 
(4.22) < G(x(p + yip + $,A c ),ip* > 



(i+u* 3 r 



d + uir a+u*r ^i)xj^ K d X 

- W) V I ^Idx + ±^x [ ^e K \ 2 dx 
(1 + ^1) / Jn 1 + Pk Jn 

+ ^T^-y ! <i>i\Ve K \ 2 dx^^A c ^lx [ e 2 K A(-A + iy^.dx 
1 + Pk Jn 2 J n 

- -^cvivly J^Ai-A + iy^idx 



VI 



1 & 

(l + u* 3 )^ 2 (1 + utf 



a^f I / e K dx 



(za&rjirjl e A K dx^j xy 2 



+ Sa^m / e K dx 



x 2 y 



arifvl / e%dx y 6 + o(3). 



Step 4. By the formula of center manifold function in the complex case in 
Theorem A.l in [5], we have 



(4.23) 

with 
(4.24) 



4>i 
$2 



*>2 



$1 



$3 



•o(3), 



$1 



= x 2 Gn + xy(Gi2 + G 2 i) + y 2 G 22 , 
(L 2 c + Ap 2 )L Xc ( = 2p 2 [(x 2 - y 2 )(G 22 - G n ) - 2a;y(G 12 + G 21 )] 



= p [(y 2 - x 2 )(G 12 + G 21 ) + 2xy(G n - G 



22 ) 



Here G,j = G 2 ( 1 I r \ A c ) with vf 1 = <^ and * 2 = -0, and G 2 is as defined in Step 
3. Namely 

Gio = 



-A c V(*iV(-A + J) _1 *i) 




+ (l+«3) 2 



,iOT l W l l| 'l 



with *f = T\e K , 1 < i, Z < 2, and 

(4.25) r}=a, r 2 = a r 2 = m , r 2 = r; 2! 
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which are given by (|4.20p . 

Direct calculation shows that 



(4.26) 



G* 



■K 



V(e K Ve K ). ., 
Pk \ 11 



2 



-arfrj 



For simplicity, we only consider the case where if 
we can see that 



(if 1,0). In this case, by flSIBJ) 
e l- = 7:( e o + e 2 if), V(e K Ve K ) = -p K &2K- 



1 



Then, by (|4~2"rJ)) . we have 
(4.27) 



eo 



<1K 



Z2K, 1 < i, j < 2, 



where 

(4.28) 

Let 
(4.29) 



•2K 



9% 



aT\T{ 



(1 T- "3 

pifA c rirj 



(1 + US)3 



- 3aulT\T{ 



1 



2A 



+ Pk 



^2 



eo 



*> 2 fe f 



e 2 A, 



Then it follows from (|4724|) and (|4T27|) that 
(4.30) 



</?22 
^32 



= -Bn 



= B. 



2K 



9°n 



2K 

9n 



xy 



xy 



/i? 2 + 

9°12 + 



52°1 

f<7 2 2 f 



1< fc < 3. 



"22 
52°2 



l2A 
"22 



2p 2 B^(B 2 



2K 



Vi)- 1 

+ 4p 2 /) 



(x 2 



ft 

"22 
ff2°2 



(* 2 -?/) 



- ft 
9n 

h 2K _ h 2K 
"22 "11 



2xy 



P(B 2 



P(B 2 



2 A' 



4P 2 /)" 1 
+ 4P 2 /)- 1 



(v 2 



h° 12 
9°i2 



^i 
<?2°i 



2A 

9n 



2xy 



2xy 



512 + 321 



9l? 



h° - h° 
"11 "22 

3?i-32°2 



2A i u2K 
"21 



hf 2 
9 2 ¥ 



2xy 



h 2K 

"li 
9 li 



h 2K 
"22 

9% 



Al c with At as defined by gl§. By (|4~2"3l> and (|4T2"9l we obtain an 



where B k — — ^ ^ fe 
explicit expression of $ as follows 

$i = (vSi 

$2 = (v»?a + <^22 + ^3 2 )eo 



(4.31) 



<P3l) e 



2A 



2A 
12 



,„2A 

,„2A 
^22 



<P 2 f ) e 2A 

2A\ 
</?32 ) e 2A 



•o(2), 
o(2). 
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Here, by (|4.30p . f|4.28[) . and (|4.25|) . are 2-order homogeneous functions of (x,y), 
with the coefficients depending explicitly on the parameters defined in (|2.4|) . 



Step 5. Finally, inserting (14.31ft into (|4.21l) and (|4.22[) . we can write (|4.18[) in 
the following form 

^ = -pty + anx 3 + a 12 x 2 y + a 13 xy 2 + a u y 3 + o(4), 

^ = px + a 2 ix 3 + a 2 2X 2 y + a 2 3xy 2 + a 24 y 3 + o(4). 
Then, by Theorem 2.4.5 in [7], the parameter 6o in Theorem 14. 2 1 is obtained by 

6 = 3an + 3a 2 4 + ai2 + 023, 
where an, a 2 4, ai2, ^23 can be explicitly expressed in the terms in fl4.20|) - (14 ,22|) . 



4.5. Transition for the system (|2.5[) . We are now in a position to discuss the 
transition of (|2 . 5|) . With the translation (14.21) . the system (|2.5[) is rewritten in the 
following form 

dui „ 2 * a aiiitt + 3 , . 

— = nAui - 2au 1 ut - u l Au 2 + ^ - ^ 2 + g{u), 



— - = Au 2 -u 2 + Ami, 
at 



du 
On 



= 0, 



u(0) = uq, 

where g(u) is as in (|4.4[) . Here the notation u stands for three-component unknown: 

u = (ui,u 2l u 3 ). 

Let 

' pA - 2auf -utA IT ^| F 
L\u = I A A - 1 0' 



-5ul rA - 5u{ 







) 


U2 




V «3 / 



Then, all eigenvalues /3^(A) and eigenvectors of La satisfy 

f C: \ ( & 

V £fc3 ' \ £fc3 

with 

^fc = (tii e k,£,i 2 e k ,e k3 ek), 
and efe as in (|3.6[) . Z?^ is a 3 x 3 matrix given by 

-(pp k + 2auf) u\p k (\Tu\)'i 
D k=\ A -(Pk + 1) 



-<y«s -(r Pfe + 0Ul ) 
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We introduce the following three parameters: 

—tiD k = fip k + 2au* 2 + pk + 1 + rp k + Sul, 



Bl = det 



det 



Ct 



-{pp k + 2au\ 2 ) u\p k 

A -(p k + 1) 

(pk + l)[rpk + Sul), 

detD^ = (p Pk + 2auf)(p k + l)(rp k + Sul) 



~(pp k + 2auf) 
-Su% 



-(rpk + Sul) 



-ulp k \{rpk + Su{) 



(i + u * 3 y 



;5ul(p k 



1). 



By the Routh-Hurwitz theorem, we know that all eigenvalues (3 3 k of Dl have negative 
real parts if and only if 



(4.33) A x k > 0, A x k B x - C x > 0, 

Let A c and K = (Ki,K 2 ) satisfy 



Ct > 0. 



(4.34) 



< c > 0, 
A K k > 0, 



AtBt 



C 



Ac 



0. 



c$? > o, 



A" 

c£° > o, c Ka > o, vfc with Pk ^ p K . 



Then A r satisfies that 



1 



(4.35) A c = inf [(/i + r)p k + 2auf + Sul 

p k p k u x 

x (r + l)p k + Sul + l + 



a5 



{p P k + 2au\ 2 + p k + 1)(1 + utY 

and pk satisfies (|4.35[) . In particular, under the condition (|4.34[) , there is a pair of 
complex eigenvalues (3 K (\) and (3 K (X) of D K , such that 



(4.36) 



< 0, A < A c , 
Ref3 K '(\){ =0, A = A C , 



> 0, A > A r 



and the other eigenvalues j3 k (\) of L\ satisfy 



(4.37) 



Re Pi(A c ) < 0, Vfc with p k ^ pk, and 1 < j < 3, 
P 3 K (A C ) < 0. 



Let A c and K* = (K{,K^) satisfy 



A%. > 0, 



(4.38) 

Then A c is given by 
(4.39) X r _ = inf 



C K \ > 0, C x % = 0, 



A k c > 0, A k c B k c - Cl° > 0, C x " > 0, Vfc with p k f pk- ■ 



(pk + 1) 



Pk p k u x 



ppk + 2aul 



*2 



a S 



(l+u* 3 ) 2 (rp k + Sul) 



24 



MA AND WANG 



and A c arrives its minimal at pa*- From the Routh-Hurwitz criterion (|4.33|) . we 
deduce that with (|4.38p there is a real eigenvalue (3^, (A) of D^, satisfies 



4 f Re/3^(A C )<0, j — 2, 3, 

' \ Ref3 J k (\ c ) < 0, Vfc G Z 2 with p fc ^ px* and 1 < j < 3. 

It is clear that (|4~36l) and (|4~37| hold true as A c < A c , and (|439)) - (|4~40)) hold 



true as A c < A c . Hence, we have the following transition theorem for (|4.32l) . 

Theorem 4.3. Let A c and X c be given by ^4-35\ l and respectively. Then, 

Assertions (1) and (2) of Theorem \4-2\ hold true for the system \^.32^ . 



5.1. Biological significance of transition theorems. Pattern formation is one 
of the characteristics for bacteria chemotaxis, and is fully characterized by the 
dynamic transitions. Theorems 13.11 - ^751 tell us that the nondimensional parameter 
A, given by 

(5.1) X= V^X 7 

plays a crucial role to determine the dynamic transition and pattern formation. 
Actually, the key factor in (|5.ip is the product of the chemotactic coefficient \ an d 
the production rate r\ : x r i > which depends on the type of bacteria. When A is less 
than some critical value A c , the uniform distribution of biological individuals is a 
stable state. When A exceeds A c , the bacteria cells aggregate to form more complex 
and stable patterns. 

As seen in (|3.11[) . (|4.10l) . (14.11[) and (14.35[) . under different biological conditions, 
the critical parameter A c takes different forms and values. But, a general formula 
for A c is of the following type: 



where pk are taken as the eigenvalues of — A with the Neumann boundary condition. 
When Q is a rectangular region, pk are given by (|3.6I) . and the coefficients dj 
(1 < j < 4), bo, b\ > depend on the parameters in (|2.4I) . with 

ao, Oi, a 2 , b , h > 0, a 3 ,a 4 >0. 

In particular, for the system with rich nutrient supplies, (|5.2p becomes 



(4.40) 




5. Biological Conclusions 



(5.2) 




(5.3) 



A c = ao + inf a x p k H . 

Pk [ Pk. 

The eigenvalues pk, depending on the geometry of f2, satisfy 
on / = po < pi < ■ ■ ■ < pk < ■ ■ ■ , Pk -> oo as k ->• oo, 




where L is the length scale of il. 
We infer from jO]) and jOJ) that 



A c -> oo as |0| ->• (L^O). 
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It implies that when the container fi is small, the homogenous state is state and 
there is no pattern formation of bacteria under any biological conditions. 



5.2. Spatiotemporal oscillation. Theorems 14.21 and 14.31 show that there are two 
critical parameters A c and A c , such that if A c < A c , the patterns formed by biologi- 
cal organisms are steady, as exhibited by many experimental results, and if A c < A c 
a spatial-temporal oscillatory behavior takes place. 

For the case with rich nutrient, 

u 1 = 1, it 3 = oo. 

In this situation, A c in (14.111) is reduced to (|3.8[) . and obviously we have that 

A c < A c for both (|4~T0f and (jOSjl . 

and the dynamic transition and pattern formation are determined by Theorems l3.ll 
and 13.21 Hence there is no spatiotemporal oscillations for the rich nutrient case, 
and the time periodic oscillation of chemotaxis occurs only for the case where the 
nutrient is moderately supplied. 
In particular, if /i, r = 0, and 

8 2 u* 2 {\ + u* 3 ) 2 < aS , 

then for A c defined by (|4.10j) and (|4.35|) . we have 

A c < A c . 

In this spatial-temporal oscillation pattern are expected for A > A c . 

5.3. Transition types. One of the most important aspects of the study for phase 
transitions is to determine the transition types for a given system. The main theo- 
rems in this article provide precise information on the transition types. In all cases, 
types are precisely determined by the sign of some non dimensional parameters; see 
b, bo andfei respectively in the main theorems. Hence a global phase diagram can 
be obtained easily by setting the related parameter to be zero. 

For example, when fl = (0,Li) is one-dimensional or when K = (Ki,0) (resp. 
K = (0, K2)), the parameter b in (|3 . 14[) can be simplified into the following form 



(5.4) 6=2 



—ZppK + 9a — 



(2fipK + a){2p,\ 2 K + 28apK + 4a — p,p K ) 



{(J-P2K + 2a){p 2 K + 1) - P2K Ac 

For a non-growth system, a = 0,K = (1,0), A c = /j,(pk+1). Then, (|5.4[) becomes 

2 

(5.5) 6=|(1-20A 1 ), Ai = ^j, 

and A = with a — jij j n u\dx. It follows from (|5.5I) that 

< if Li < 2v / 5tt, 
> if Li > 2V5tt. 



(5.6) b 



By Theorems 13.11 and 13.21 the phase transition of f|3 . 3[) and (|3.ip from (u, A) = 
(u* , A c ) is continuous if the length scale L\ of O is less than 2\/57r, and jump if L\ 
is bigger than 2\/5tt. 

In addition, when we take 

x{u) = V + Ztf 
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as the chemotaxis function, by Remark |4.1[ the parameter b of (|5.5|) is replaced by 

u ( 20vr 2 \ 4ku 2 tt 2 
°i = 77 1 



3 V~ L\ J L\ 

with 

" = /dW andA c = M g + l 

The above conclusion amounts to saying that for a non-growth system, the pa- 
rameter 

A = ^ a, with a = — — ■ / u\dx, 

T2K2 |"| 

is proportional to the average density a of initial condition of u\(u\ is conservation). 

Hence, the biological individual is in a homogenous distribution state provided 



tpdx < -^fJ- ( 72 + 1 ) , = ui(0), 



1 / ^ / r2fc2 

and the bacteria will aggregate to form numbers of high density regions provided 

1 f , r 2 k 2 ( 7T 2 
(5.7) (pd x> -7— -/i to + 1 



Moreover, under the condition (|5.7|) . if the scale Li of is smaller than some critical 
value L c (in (|5.6p L c = 2\/57r), i.e. ii < L c , the continuous transition implies that 
there is only one high density region of bacteria to be formed, and if L\ > L c then 
the jump transition expects a large number of high density regions to appear. 

5.4. Pattern formation. As mentioned before, the pattern formation behavior is 
dictated by the dynamic transition of the system. In this article, we studied the 
formation of two type patterns-the lamella and the rectangular patterns, although 
the approach can be generalized to study the formation of other more complex 
patterns. 

For a growth system, the critical parameter A c takes its value at some eigenvalue 
p K of -A for K = (Ki,K 2 ), as shown by (13.111) and (I4.11[) . From the pattern 
formation point of view, for the Type-I transition, the patterns described by the 
transition solutions in thee main theorems are either lamella or rectangular: 

lamella pattern for K\K 2 = 0, 

rectangular pattern for K\K% 7^ 0. 

In the case where b > 0, the system undergoes a more drastic change. As A* < 
A < A c , the homogeneous state, the new patterns v 2 and are metastable. For 
A > A c , the system undergoes transitions to more complex patterns away from the 
basic homogeneous state. 
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